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ABSTRACT 

Constraining the properties of Population III (Pop III) stars will be very challenging because 
they reside in small galaxies at high redshift which will be difficult to directly detect. In this 
paper, we suggest that intensity mapping may be a promising method to study Pop III stars. 
Intensity mapping is a technique proposed to measure large-scale fluctuations of galaxy line 
emission in three dimensions without resolving individual sources. This technique is well 
suited for observing many faint galaxies because it can measure their cumulative emission 
even if they cannot be directly detected. We focus on intensity mapping of He II recombina¬ 
tion lines. These lines are much stronger in Pop III stars than Pop II stars because the harder 
spectra of Pop III stars are expected to produce many He II ionizing photons. Measuring 
the Hell 1640 A intensity mapping signal, along with the signals from other lines such as 
Lya, Ha, and metal lines, could give constraints on the initial mass function (IMF) and star 
formation rate density (SFRD) of Pop III stars as a function of redshift. To demonstrate the 
feasibility of these observations, we estimate the strength of the Pop III He II 1640 A intensity 
mapping signal from z = 10 — 20. We show that at z ~ 10, the signal could be measured 
accurately by two different hypothetical future instruments, one which cross-correlates He II 
1640 A with CO(l-O) line emission from galaxies and the other with 21 cm emission from the 
intergalactic medium (IGM). 
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1 INTRODUCTION 


There has been much theoretical work examining Pop III stars 
assembled from primordial gas in the early Universe. The first 
Pop III stars are expected to form in (10 5 — 10 s )A fp , miniha- 
los that undergo molecular hydrogen c ooling dHaiman et al .11 1 9961 : 
iTegmark et alii 19971 : 1 Abel et al] 120021. As the Universe evolves, 
stars produce a cosmological background of Lyman-Werner 
(LW) radiation which dissociates molecular hydrogen, eventu- 
ally preventing addi ti onal star formation in low-mass minihalos 
(HaimanetaL 19971: Machacek et al. 2001Wise&Abel 20071 


O’Shea & Normanl 120081: IWolcott-Green et al.11201 ll: Ivisbal et al] 


20141) . Star formation is then restricted to larger halos, with masses 


near that of atomic cooling halos (i.e. halos with virial temperature 
above T v i r ss 10 4 K). Three-dimensional hydrodynamical simu¬ 
lations have been utilized to understand the formation of Po£_III 
stars both in minihalos and atomic cooling halos (e.g. Stacy et al. 

J I2012I: 


,20101: Clark et al 


Wise et al. 120121: 


2011; Hosokawaet al. 2012: iGreif et alJ 


Hirano et al. 20141) . Although these simulations 


cannot make predictions for the properties of Pop III stars, they 
have provided estimates for the sizes of the gas clumps forming 
as a result of fragmentation. The results suggest that the IMF of 
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Pop III stars may have been top-heavy, but the IMF remains highly 
uncertain. 


Unfortunately, it will be extremely challenging to test predic¬ 
tions for Pop III stars by observing them directly. This is because 
they are predicted to form in very small halos at great distances. 
In fact, even the James Webb Space Telescope (JWST) llGardnetl 
l200d) may not be able to detect Pop III stars in atomic cooling ha¬ 
los at the epoch of reionziation (z ~ 10) dZackrisson et al.lf201 li) . 
Pop III stars in minihalos at earlier times will be even more dif¬ 
ficult to observe. While it will be very interesting to study larger 
systems at high redshift, these galaxies are likely to have been en¬ 
riched with metals from previous generations of star formation and 
therefore may not host significant Pop III star formation. 


In this paper, we propose a new technique to constrain the 
properties of Pop III stars, intensity mapping of helium recom¬ 
bination lines. Intensity mapping is a technique to measure the 
large-scale three-dimensional clustering of galaxies by observing 
emission line fl uctuations as a function of wavelength and posi¬ 
tion on the sky d Visbal & Loebll2010l : Ivisbal et alj|201 1| . see also 
iRighi et all 120081). and has b e en proposed for a number of dif- 
ferent line s dGong et ai]|201 ll: lOdz et al.ll201 ll: IGong et alJlioTl 


ISilva et alJ 120131: IGong et alj|2013l: IPullen et alJl2014h . Intensity 
mapping does not require resolving individual sources, instead it 
measures the cumulative line emission from all galaxies in coarse 
three-dimensional pixels. This makes intensity mapping a promis- 
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2 E. Visbal et al. 


ing technique to study the small systems that host Pop III star 
formation. Even though individual Pop III galaxies may be unde¬ 
tectable, their large numbers will result in a cumulative signal that 
could be measured through intensity mapping. 

We focus on He II lines (and in particular He II 1640 A) be¬ 
cause, compared to metal-enriched (“Pop II”) stars, the hard spectra 
of Pop III stars are expected to produce many more He II ioniz¬ 
ing photons relative to HI ionizing photons, leading to strong re¬ 
combination lines ( Schaereill2002|j2003t) . For this reason, Pop III 
stars are likely to dominate the He II intensity mapping signal (see 
Section 5). Studying Pop III stars with He II lines has bee n con ¬ 
sidered in the context of individual sources |Tumlinson & Shulll 


2000; lOh et al .1 [201)71: ICai et al J[2oTTL I20l4 ICassata et al.ll2013l) 


but as mentioned above, it will be very difficult to detect the 
smallest Pop III sources, even with future telescopes. He II lines 


down to a ~ 4 (e.g. 

Dawson et al. 2004; Jimenez & Haiman 2006; 

Eldridge & Stanwav 

20 lj, with current results suggesting that the 


contribution of Pop III stars to the total SFRD is at most at the per¬ 
cent level feheng et al.ll2013l) . Measuring the He II intensity map¬ 
ping signal, along with the intensity mapping signal of other lines 
such as Lya, Ha, and metal lines, will constrain the amount of Pop 
III star formation as a function of redshift and the hardness of their 
composite spectra which depends on the IMF. 

This paper is structured as follows. In Section 2, we give a 
brief review of the intensity mapping technique. We discuss how 
the intensity mapping signal can be used to constrain the properties 
of Pop III stars in Section 3. In Sections 4, 5, and 6, we estimate 
the strength of the He II 1640 A intensity mapping signal from Pop 
III stars, Pop II stars, and quasars, respectively. In Section 7, we 
describe two hypothetical future intensity mapping experiments, 
one which cross-correlates Hell 1640 A with CO(l-O) (2610 /im) 
emission and the other with 21 cm emission from the IGM. As we 
discuss below, cross-correlation is necessary to remove contamina¬ 
tion from other emission lines originating from galaxies at different 
redshifts. Finally, we discuss our results and conclusions in Sec¬ 
tion 8. Throughout we assume a ACDM cosmology consistent with 
the latest constraints from Planck iAde et al]|20 14h : Ha = 0.68, 
Sim = 0.32, Hb = 0.049, h = 0.67, ag = 0.83, and n s = 0.96. 


2 INTENSITY MAPPING REVIEW 

Here we review technical details of the i ntensity mappin g tech¬ 
nique. For a more thorough discussion see I Visbal & Loebl <20101) : 
Visbal et al.l d20 1 ill . Intensity mapping requires measuring fluctua¬ 
tions in line intensity as a function of position on the sky and fre¬ 
quency, corresponding to a three-dimensional position in redshift + 
sky position space. The total intensity mapping fluctuation signal 
for a particular line, after removing the spectrally smooth compo¬ 
nents including galaxy continua and foreground/background radia¬ 
tion, is given by 

AS to t(r) = AS(r) + ^ A5' CO nt,i(r') + ASWoise, (1) 

i 

where A S is the signal from the target line at the target redshift 
+ sky position pixel r, the AS CO nt,i’s are the contaminating emis¬ 
sion lines appearing at the target frequency from unrelated galaxies 
at different redshifts, and ASNoise is the detector noise. Assuming 
galaxies are found in dark matter halos which trace the large-scale 
cosmological density field, spatial fluctuations in galaxy line inten¬ 
sity are given by 


AS = SbS + ASpoisa, (2) 

where S is the mean line signal (proportional to the total line emis¬ 
sion per volume), b is the mean line luminosity-weighted galaxy 
bias, <5 is the cosmological dark matter overdensity, and ASWss is 
the variation due to Poisson fluctuations in the number of galaxies 
at fixed 5. Throughout we approximate the mean bias as 

fdMb(M)L(M)& 

JdML(M ’ 

where L(M ) is the line luminosity from halos of mass M, b(M) 
is the Sheth-Tormen halo bias, and -fe is the halo mass function 

dSheth & Tormenll 19991). _ _ 

As discussed in Visbal & Loebl I2010f) . it is possible to isolate 
emission from a particular target redshift by cross-correlating dif¬ 
ferent lines from the same volume. The cross-correlation function 
is defined by 

£i, 2 (r) = (A5 , i(x)A5 , 2 (x + r)) (4) 


where the subscripts denote different emission lines. The cross 
power spectrum is defined as the Fourier transform of the cross¬ 
correlation function 

A, 2 (k) = J d 3 r£ 1 , 2 (r)e* k ' r , (5) 

and (ignoring redshift space distortions) is given by 

A, 2 (k) = 5ifei5 2 6 2 P(k) + P Poiss , (6) 

where P( k) is the matter power spectrum and -Ppoiss is the (con¬ 
stant) contribution from the Poisson fluctuations. The error on a 
measurement of the power spectrum for one k-rnode is given by 

^Pi,2 = 2 (-ft,2 + Ptot,lPtot,2) , ( 7 ) 

where Ptot,i and Ptot,2 are the total power s pectr a for each d ata set 
being cross-correlated (see Appendix A of I Visbal & Loebl (201(1) 
for a derivation of this formula). This is given by 

Ptot.l = Pi + PNoise + E -fcont ,i 5 (8) 


where Pi is the power spectrum of the target line, Pn 0 ; se is the 
contribution from detector noise, and the P CO nt,i’s are the power 
spectra from contaminating lines at different redshifts. These lines 
come from galaxies whose positions are assumed to be uncorre¬ 
lated with those of the target galaxies because they are at different 
redshifts separated by very large distances (see below). 


3 INTERPRETING THE INTENSITY MAPPING SIGNAL 

In order to extract the most information from intensity mapping 
experiments it will be useful to measure the signal from a number 
of different lines. By taking the ratios of different lines’ clustering 
cross power spectra (i.e. the first term in Eqn. [6jl, one can find Sb 
for each line as a function of redshift. Additionally, the mean signal 
and bias could potentially be separated by examining the angular 
dependence of the power spectrum due to redshift space distortions. 

To study Pop III stars it will be useful to measure the inten¬ 
sity mapping signal of He II 1640 A, Lya, Ha, and a number of 
different metal lines. As discussed below, two possible (but most 
likely sub-dominant) non-Pop III contributions to the He II signal 
come from Wolf-Rayet (WR) stars and (mini-) quasars. By measur¬ 
ing the intensity mapping signal of metal lines expected from WR 
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stars or quasars in addition to He II, it will be possible to verify 
that Pop III stars dominate the He II signal and potentially subtract 
other contributions. It may also be possible to use metal lines to 
estimate the total Ha or Lya emission from Pop II stars and sub¬ 
tract it away, isolating the Pop III signal. One could then infer the 
total He II 1640 A and Ha signal from Pop III stars as a function 
of redshift. The ratio of these lines is a measure of the hardness 
of Pop III s pectra which de pends on the IMF (see lOh et alj|200ll 
and fig. 5 in Schaerer 2003). The Lya line could be combined in 
a similar way, but would yield more model dependent results due 
to radiative transfer effects in the IGM tha t strongly modi f y the ap¬ 
parent Lya lineshape (e.e. lHaima ^ <20021) ; ISantos|j2004l) JPiikstral 
<2014l) ). The He II signal as a function of redshift would give infor¬ 
mation about how the SFRD of Pop III stars evolves over cosmic 
time. 


4 HE II 1640 A LINE SIGNAL FROM POP III STARS 


In this section, we estimate the strength of the signal from Pop III 
stars to demonstrate the feasibility of He II 1640 A intensity map¬ 
ping. We assume that the line luminosity is given by 


LHell = R X SFR X (1 — /esc) , 


(9) 


where SFR is the star formation rate, R is a constant that depends 
on the IMF (and the metallicity if the stars are not metal-free), 
and / e sc is the escape fraction of Hell-ionizing (>54.4eV) photons, 
which we assume to be much less than unity. For constant star for¬ 
mation, zero metallicity, and a Salpeter IMF with a mass range of 
50 - 500, 1 - 500, or 1 - 100 M 0 , R has a value of 1. 6 x 10 8 
2.5 x 10 7 , or 4.6 x 10 6 Lq/(M@ yr _1 ), respectively dSchaererl 
120031) . 


In the examples discussed below, we consider He II inten¬ 
sity mapping at redshifts z > 10. We note that this is a con¬ 
servative choice, as Pop III stars could ex i st at redshifts as low 
as z ~ 3 — 4 1 Jimenez & Haimaiil 120061 : IWvithe & Cenl 1 20071 : 
iFumagalli et al.ll201ll) potentially producing a stronger signal, ex¬ 
tending down to lower observed wavelengths. We assume that the 
LW background is strong enough to suppress star formation in 
minihalos and focus on atomic cooling halos (but will address a 
possible extra contribution from minihalos below). For dark matter 
halos above the atomic cooling threshold we approximate the SFR 


SFR = 


/* 


fib M 

flm Cdutyltt^) 


( 10 ) 


where M is halo mass, /» is the star formation efficiency, eduty is 
the duty cycle, and t-n(z) is the Hubble time. To compute the total 
signal, we assume that Pop III stars are formed in halos with M = 
(1 — 2)M coo i, where M coo 1 is the atomic cooling mass (M coo 1 = 
3 x 10 7 ((1 + z)/ 11) 1 ' 5 M 0 corresponding to TLir = 10 4 K for 
neutral gas). Much larger halos are likely to have been enriched 
with metals from stars formed in their progenitor halos. Given these 
assumptions, the mean Pop III He II signal (i.e. the specific surface 
brightness, in units of ergs s _1 cm -2 Hz -1 sr _1 and proportional 
to line luminosity per volume of space) is 


/ 
J N 


o _ ( 2M “ o1 dn ~ r >2 

Snell - / dM D ,j £ duty dM vD A , 


(ID 


-^cool 


where Dl is the cosmological luminosity distance, Da is the angu¬ 
lar diameter distance, and y = is the corresponding comoving 
distance (x) per change in observed frequency. Note that the mean 



Figure 1. Mean Hell 1640 A specific surface brightness versus redshift 
(corresponding to an observed wavelength of 1640 A x (1 + z )) for differ¬ 
ent Pop III IMFs and /* =0.1. The IMFs are Salpeter with mass ranges of 
(1 — 1 OO)M 0 (bottom solid curve), (1—5OO)M0 (middle solid curve), and 
(50 — 500) Mq (top solid curve). The dotted curve is the signal from mini¬ 
halos assuming no LW background (assuming a Salpeter IMF with mass 
range (50 — 500) Mq). The dot-dashed curve is the upper limit from (Pop 
II) WR stars discussed in Section 5 and the dashed curve is the contribution 
from quasars estimated in Section 6. 


signal does not depend on eduty. We plot the signal as a function 
of redshift for several different IMFs and /* =0.1 in FigureQ] For 
comparison, we also plot the signal for Pop III stars in minihalos 
assuming no LW background. We compute this using Eqn. El but 
replace the atomic cooling mass with the cooling mass for molec¬ 
ular cooling 2 _MH,= 2.5 x 10 5 3 ^ 2 Mp) 1 Machacek et all 

1200 it IWise & Abelll2007l: Id’ Shea & Normadl2008h . Note that not 


including LW feedback results in a large overestimate below the 
redshift where this effect becomes important. In Figure [2] we plot 
the SFRD of Pop III stars corresponding to the He II specific sur¬ 
face brightness in FigureQ] We note that for Pop III stars in atomic 
cooling halos, the SFRD is within roughly an order of magnitude of 
that computed in lTrenti & Stiavellil 420091) when they assume mul¬ 
tiple Pop III stars can form in each dark matter halo. 

By using a constant duty cycle, combined with a snapshot 
of the halo mass function (rather than using the halo formation 
rate) Eqns. 1101 and El give a rough approximation of Simi and 
the corresponding Pop III SFRD. To test this approximation, we 
estimate the SFRD (which is proportional to Simi) with a differ¬ 
ent approach and compare the results. Numerical simulations show 
that at z ~ 10, the rate at which halos cross the atomic cooling 
threshold per volume is dn J° o1 ~ 4 Mpc -3 Jvisbal et alJl2014l) . 
Assuming /* =0.1 of the gas in these halos forms Pop III stars 
over a time period £sf and that dn J° o1 does not evolve rapidly over 
this time scale, the Pop III star formation rate density is given by 
SFRD = dn ™ Ql (/* ^-M coo i^. The SFRD computed this way 
is ~ 3.5 times greater than that implied by Eqns. Qo] and El (this 
increase is attributable to the rapid halo formation rate). Thus, the 
actual fraction of gas that is required to form stars to reach a partic¬ 
ular SHeii may be less than /* in Eqn.1101 

To check that our SFR prescription and fiducial choice of /* = 
0.1 are reasonable, we compare the amount of ionizing photons 
produced at z ~ 10 to that required to maintain reionization. The 
recombination rate per volume in a completely ionized region is 
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Figure 2. The SFRD corresponding to the He II specific surface brightness 
plotted in Figure [I] for atomic cooling halos (solid curve) and minihalos 
(dotted curve). Note that we have ignored the negative feedback on miniha¬ 
los due to LW radiation. 


B — obWigmCi 0-2) 

where qb = 2.6 x lCP 13 cm 3 s -1 is the case B recombination co¬ 
efficient, niGM is the mean hydrogen number density of the inter- 
galactic medium (IGM), and C = (niGM,Hii)/( n iGM,Hii ) 2 is the 
clumping factor of the ionized gas in the IGM. The rate of ioniza¬ 
tions per volume from star formation is given by SFRD x Q(H) x 
/esc, hi, where Q(H) is the number of ionizing photons per SFR 
and /esc ,hi is the escape fraction of HI ionizing photons. We com¬ 
pute the SFRD using Eqn. [To] and integrating the halo mass func¬ 
tion. We assume that Pop III stars with a Salpeter IMF and mass 
range of (50-500 )Mq form in halos with M = (1 — 2)M CO oi and 
that Pop II stars with a metallicity of Z = 0.0004 and a Salpeter 
IMF with mass range ( l-lOO)Mro form in larger ha los. Using the 
values of Q(H) from ISchaererl d2003l) . C = 3 dFinlator et al.l 
|20I2|) . /* = 0.1, and /esc ,hi = 0.03, we find that the rate of re¬ 
combination equals the rate of ionization if the IGM were ionized. 
Thus, our fiducial choice of /* =0.1 appears reasonable. 


5 POP II WOLF-RAYET SIGNAL 

The most significant source of He II recombination lines from Pop 
II stars is likely to come from WR stars. Here we show that this 
component is unlikely to dominate the signal at high-redshift. It is 
difficult to estimate the WR signal from populati on synthesi s mod¬ 
els, so instead we utilize observational results. ICassata et ail d2013l) 
compiled a sample of high-redshift galaxies (z = 2 — 4.6) with 
broad Hell lines consistent with those expected from WR stars. 
For each galaxy, they compare the SFR estimated from the broad¬ 
band SED to the He II 1640 A luminosity and find the best fit linear 
relation 

T/Heii,wR = 3.5 x 10 39 I —- -r ) ergs s 1 . (13) 

V M © yr- 1 / 

This can be viewed as an upper limit on the WR contribution since 
most galaxies’ spectra do not show strong He II spectral features. 

This relation is consistent with the empirical estimates of 
Q(HeII)/Q(H) (where Q s are ionization rates) inferred from 
the He II 4686A to H /3 line ratios in extragalactic HII regions 


discussed in Schaered d2003l) (see their sectio n 6.3, see also 
iGuseva eta i] |2000h . We use the stellar models of ISchaeretl d 20031) 
to translate these estimates into He II 1640 A luminosity per SFR. 
Assuming a Salpeter IMF with a mass range of M = (1 — 100) Mg 
and a constant SFR, we find that for metallicity Z = 10 -3 , 


F'Hell.WR ~ (1.4 — 5.5) x 10 3 


SFR 
Mg yr~ 


ergs s 


(14) 


The observations indicate that at higher metallicity the WR He II 
signal is much lower. 

Assuming that Pop II stars are found in dark matter halos 
above 2M coo i and using Eqn. |T3] and Eqn. [TO] we plot the up¬ 
per limit of the WR contribution to the Hell 1640 A intensity 
mapping signal in Figure Q] The Pop III component dominates at 
all redshifts for top heavy IMFs (i.e. Saltpeter with mass limits 
of (50 — 500)Mg or (1 — 500 )Mq); and it remains dominant 
at redshifts well above z > 10 for a usual Salpeter IMF (with 
(1 - 100)M©). 


6 QUASAR SIGNAL 


The spectra of quasars are harder than stellar sources and thus can 
also efficiently ionize He II. Here we estimate their contribution to 
the Hell 1640 A intensity m apping sig nal. To calculate the line 
strength, we use eqn. 2 fro m iQh et al.l ( 20011 ) . which assumes a 
spectral template from IeI vis et al.l 1 199-41) normalized such that the 
total luminosity corresponds to the Eddington limit. This gives an 
He II 1640 A luminosity of 


T/Hell, quasar = 10 41 ( ) er S S S S (15) 

where Mbh is the mass of the black hole powering the quasar (or 
mini-quasar). We assume that black hole mass relates to halo mass 
as 


Mbh = 3 x 10 


M 


5 x 10 12 Mg 


5/3 


1 +z 


Me,. 


(16) 


and compute the mean intensity mapping signal from the halo mass 
function. We plot the average signal versus redshift in FigureQ]as- 
suming that the atomic cooling mass is the minimum halo mass 
hosting quasars and that they have a duty cycle of 0.05 (a high duty 
cycle compared to low-redshift quasars is possible due to the in¬ 
creased frequency of major mergers at high redshift. fTanakall20 14l ). 
For this choice of duty cycle, Eqn. M implies that the density 
of active quasars hosting ~ 3 x 10 9 A/q or larger bl ack h oles is 
~ 10~ 9 Mpc -3 , which is consistent with observations lFanll2006l) . 
The mass and redshift scaling are expected for feedback-limited 
accretion ( iLoeb & Furlanettol2013l) . We find that the Pop III signal 
dominates for Saltpeter IMFs with mass limits of (50 — 500) Mq 
or (1 — 500) Mq at z > 10. 


7 FUTURE OBSERVATIONS 

In this section, we discuss the capability of a hypothetical futur¬ 
istic instrument to measure the Hell 1640 A intensity mapping 
signal. We consider a space-based instrument with a 2 m dish and 
background-limited sensitivity. The background radiation at the rel¬ 
evant wavelengths _is dominated by Zodiacal light and faint stars. 
We use fig. 1 of lLeinert et al .1 dl998l) to determine its strength. We 
assume that the instrument can simultaneously measure spectra of 
100 separate 4 arcmin 2 pixels with spectral resolution v /Az/ = 
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1000. This design represents an instrument which is not feasible 
with current technology, but could potentially be built in the near 
future. 

We consider two separate examples, one cross-correlating 
Hell 1640 A with CO(l-O) emission from galaxies and the other 
cross-correlating with 21 cm emission from the IGM. For both ex¬ 
amples, we compute the total signal-to-noise of the cross power 
spectrum defined as 


S_ 

N 


\ 


E 

k— modes 


P T2 


(17) 


where the sum is over the available k-modes, determined by the di¬ 
mensions and spatial resolution of the survey. We note that modes 
at the redshifts and scales we consider are in the linear regime for 
the target galaxies and thus are statistically independent. Ignoring 
redshift-space distortions (described below) and assuming the value 
of the matter power spectrum can be estimated accurately from con¬ 
straints on cosmological parameters, Jj- corresponds to the uncer¬ 
tainty on S 161 S 2&2 measured with an inverse variance-weighted 
average. We note that for the observations described below, the 
Poisson component of the signal is negligible compared to the clus¬ 
tering component for the scales we consider. 

A large source of noise in the He II cross power spectrum 
comes from emission in numerous contaminating lines (e.g. Ha) 
from foreground and background galaxies (the third term on the 
right hand side of Eqn. [ 8 j. To estimate this noise, we assume that 
galaxy line luminosity is proportional to SFR and calibrate the 
strength of the emission b ased o n the typi cal st arbur st gala xy NGC 
7714. We use Table 3A in Gonzalez-Delgado et alj i 19951) and Ta¬ 
ble 5 in iRiffel et ahl d2006l ) to determine the relative strength of 
~ 50 lines with rest frame wavelengths between 3727A and 2.2/itn. 
Contaminating galaxies are assumed to be hosted by dark matter 
halos with mass greater than 2M coo i at 2 > 6 and Eqn.[T 0 ]is used 
to determine their SFR. For lower redshift, we assume that radia¬ 
tive feedback on the IGM from reionization increases the minimum 
virial temperature of halos hosting galaxies to T v n = 10 5 K. We 
assign a SFR to these galaxies that is proportional to halo mass, 
with a normalization that sets the total SFRD (determined by inte¬ 
grating the halo mass function) to 


U + bZ (h 


(18) 


where a = 0.0118, b = 0.08, c = 3.3, and d = 5.2, which agrees 
with observations (Hopkins & Beacom 20o3). 

In order to increase the signal-to-noise of the cross power 
spectrum measurement, it will be useful to either mask or sub¬ 
tract away line emission from bright foreground galaxies in the 
He II intensity map. We consider two different conditions for re¬ 
moving these contaminating lines. In the more conservative case, 
we assume that lines are removed if their strengths (calibrated with 
NGC 7714 as described above) are more than 5 times the detec¬ 
tor noise. This effectively sets a maximum mass of contaminat¬ 
ing foreground halos (i.e. we impose an upper limit on the integral 
over halo mass for the contamination signal, analogous to Eqn. II 11 1 . 
In practice, it would be necessary to identify a few lines at a few 
different wavelengths from the same contaminating galaxy. We do 
not consider this procedure in detail, instead we assume galaxies 
with high signal-to-noise lines could be identified and removed. 
In the second, optimistic case, we assume that data from other 
telescopes could be used to mask or subtract line emission from 
bright foreground galaxies. In particular, we remove any galaxies 


which could be observed with the Large Synoptic Survey Tele- 
scoptfl (LSST) “deep drilling” fields. Galaxies with UV magnitude 
less than biab = 28 are removed. The UV luminosity is computed 
with the following relation dKennicuttll 1998l) . 


SFR 

[Mq yr-i] 


= 1.4 x 1CT 28 


L v 

ergs s _1 Hz -1 


(19) 


In our signal-to-noise calculations, we also include the linear 
change to the po wer spectrum from redshift-space distortions (see 
the Appendix of lGong et alj|2014l ). which effectively changes the 
bias from b —> b + cos 2 6 , where D is the growth function, 

a is the scale factor, and 9 is the angle between k and the line 
of sight. This is modified slightly for background and foreground 
galaxies due to the varying relationship between distance and cor¬ 
responding angle/observed frequency. We also note that the matter 
power spectrum can enter the non-linear regim e for low-redshift 
foreground galaxies. We use the fit from lTakahashi et~akl d 20121) to 
compute the increase in power due to non-linearity. 


7.1 CO(l-0)-He II cross-correlation 


For our first example, we consider a measurement of the CO(l-O) 
(2610 /.im)-He II 1640 A cross power spectrum. To estimate the 
strength of the CO signal, we assume that CO line emission comes 
from halos above 2M coo i (while halos below this mass host Pop III 
star formation) and the line lu minosity is proport ional to halo mass. 
We use the normalization from lLidz et al.l J201 ll ). 

^co(i-o) = 2.8 x 10 3 

and assume that CO emission comes from a fraction of halos cor¬ 
responding to a duty cycle eauty — 0.1. These assumptions give a 
mean signal of Sco(i-o) = 4.5 fj,K (15 Jy sr _1 ) at 2 = 10. This 
is similar to the results of other studies, but corresponds to opti¬ 
mistic estimates for a very uncertain signal. The true signal could 
be lower if, for example, low-mass galaxies near 2M coo \ do not 
have sufficient metals to produce a strong CO(l-O) line. 

We examine the potential of the Hell 1640 A intensity map¬ 
ping instrument described above and a separate CO experiment. 
We assume the CO experiment matches the 2 arcmin angular res¬ 
olution of the He II instrument and that both intensity maps have 
a total sky coverage of 400 deg 2 over the entire course of the ob¬ 
servation (combined as a mosaic from many pointings on the sky). 
We also assume matching spectral resolution of v/ Au = 1000 and 
that each observation spans a total bandwidth of B = O.lu. 

For the CO(l-O) survey, we consider an interferometer with 
detector noise 


on 


T sys 1 

%/Ai/ffl e ld /cover 


( 21 ) 


where the system temperature is taken to be T sys = 30 K and 
ffield is the integration time per field of view (assuming the ob¬ 
servation is a mosaic of many such fields). The covering factor 
/cover = N^D 2 / ax , where JV a is the number of antennae, D 
is the size of each antenna, and Umax is the maximum base line, 
which sets the angular resolution. Assuming D = 0.4 m (corre¬ 
sponding to a 4° field of view at 2 = 10), Umax = 50 m (corre¬ 
sponding to 2 arcmin resolution at 2 = 10 ), JV a = 10 , 000 , and a 
total integration time of 1 yr, we find on ~ 10 /rK (33 Jy sr _1 ). 
We use this value of ai v for all of our calculations. This represents 


1 http://www.lsst.org/ 
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Figure 3. Top panel: signal-to-noise for detecting the He II1640 A-CO(l-O) 
cross power spectrum for the hypothetical surveys described in Section 7.1 
(see Eqn. ED- We assume a Salpeter IMF with mass range of (50-500 )Mq 
(a mass range of (1 — 500) Mq would result in a signal-to-noise ~ 6 times 
lower). The solid and dashed lines are for our conservative and optimistic 
prescriptions for removing foreground galaxies. For reference, the dotted 
line is the signal-to-noise with no contaminating lines from different red- 
shifts. Bottom Panel: The contribution to the signal-to-noise squared as a 
function of scale for the conservative contamination case at z = 10. The 
contribution is the total from all modes in | k |-bins. Most of the information 
comes from scales corresponding to k ~ (0.2 — 0.7) Mpc 1 . 


a futuristic, second generation intensity mapping experiment. Note 
that this formula assumes uniform coverage in Fourier space and is 
a conservative esti mate, s omewhat fewer antennae may be needed 
for this sensitivity dLidz et al.l201 lh . The noise assumed here is in¬ 
tended as a rough estimate, future work will be required to optimize 
the design of CO intensity mapping instruments. For the He II sur¬ 
vey, we determine the noise by assuming that the entire observation 
consists of one year of integration time. 

In Figure [3] we plot the signal-to-noise of the cross power 
spectrum for this measurement as a function of redshift assuming 
a Salpeter IMF with a mass range of (50-500 )Mq. The signal- 
to-noise would be roughly 6 times lower for a mass range of 
(1 — 500 )Mq. The Pop III He II 1640 A signal could potentially 
be accurately measured at z = 10 and beyond. 


7.2 21 cm-He II cross-correlation 

Next we consider a future measurement of the He II 1640 A-21 cm 
cross power spectrum. To estimate the value of the power spectrum 
we use the publicly available code 21cmFAST dMesinger et al .1 
1201 lb . We compute the ionization field at a = 10 assuming a mini¬ 
mum halo mass of ionizing sources equal to M coo \ and an ionizing 
efficiency that sets a region as ionized if its collapsed fraction is 
greater than 1/12. This leads to a neutral fraction of xh = 0.68 at 
2 = 10. Assuming the IGM temperature is much greater than CMB 
temperature and ignoring peculiar velocities, the 21 cm fluctuation 
signal is given by 

AT 21 = 30z h (1 + 5b) mK, (22) 

where ih is the neutral fraction as a function of position and the 
baryon overdensity, 5b. is assumed to equal the total matter over¬ 
density, 5. We then compute the cross power spectrum by taking the 
discrete Fourier transform of AT 21 and the He II 1640 A intensity 
mapping signal, ASHeii = Sn e iibS. The 2 = 10 power spectrum 
is shown in Figure [4] Note that the HeII-21 cm cross power spec¬ 
trum is not of the form of Eqn. [ 6 ] However, the ratio of the cross 
power spectra can still be taken for different galaxy lines to give Sb 
as a function of redshift. For simplicity we do not consider redshift 
space distortions for the He 11-21 cm power spectrum. 

We calculate the signal-to-noise of the cross power spec¬ 
trum for the hypothetical He II instrument described above and a 
next generation 21 cm experiment. For the 21 cm experiment, we 
consider a scaled up version of the Murchison Widefield Arra)0 
( MWA ). We assume an interferometer with 5000 tiles each with ef¬ 
fective area A e = 16 m 2 . The tiles are distributed in a constant core 
up to a radius of 60 m and fall off at larger radii as n tu e oc r~ 2 . 
The layout determines the baseline distribution, which controls k- 
mode coverage in the direction perpendicular to the line of sight. 
The total 21_ cm power spectrum (including detector noise) is given 
by dMcOuinn et al. 112000) 

-fl,<iO + 0ls(%)’ i Aj, ,23) 

where P 21 is the cosmological 21 cm signal, A = (1 + 2 ) x 21 cm, 
t = 3 yr is the total observation time (1 yr total integration time 
is assumed for the Hell observation), and n(k±) is the baseline 
density. The system temperature is dominated by the sky temper¬ 
ature (mostly from Gala ctic synchrotron emission ) and given by 
T sys = 240 (ipj-) 2 ° K J Rosters & Bowmanl2008b . The total field 
of view of the observation is taken to be 1000 deg 2 and we assume 
that the angular resolution of the He II instrument is 3 arcmin in 
this case to better match the 21 cm experiment. The total signal-to- 
noise at 2 = 10. given by equationl771is S/N = 21, assuming the 
optimistic condition for removing foreground galaxies in the He II 
map. We plot the contribution from different scales in Figure [5] 


8 CONCLUSIONS 

In this paper, we argue that intensity mapping is a promising tech¬ 
nique to constrain the properties of Pop III stars because it can 
measure the cumulative line emission from many sources that are 
too faint to be resolved individually. This is particularly useful for 
studying Pop III stars because they are expected to be found in very 

2 http://www.mwatelescope.org/ 
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Figure 4. The absolute value of the He II 1640 A-21 cm cross power spec¬ 
trum at z = 10. The sign of the cross power spectrum is negative for 
k < 0.7 Mpc~ 3 , where the dip occurs (which is related to the typical 
size of ionized bubbles in the IGM). For Pop III stars, we have assumed a 
Salpeter IMF with a mass range of (50 — 500)Mq. 



Figure 5. The contribution to the signal-to-noise squared as a function of 
scale for the HeII-21 cm cross power spectrum at z = 10. The total signal- 
to-noise is -S. = 21. For Pop III stars we have assumed a Salpeter IMF 
with a mass range of (50 — 500) Mg. We use the optimistic consideration 
described in the text for contaminating line removal in the He II intensity 
map. 

faint galaxies at high redshift which will be difficult to observe in¬ 
dividually even with future telescopes such as JWST. 

We focus on intensity mapping of He II recombination lines 
(the 1640 A line in particular) because the hard spectrum pre¬ 
dicted for Pop III stars makes it likely that their contribution will be 
dominant compared to Pop II stars or quasars (see Sections 5 and 
6 ). Measuring the Hell 1640 A intensity map ping signal a s well 
the signal from other lines such as Ha, Lya dPullen et alj|2014l : 
iSilva et alj[2013l) . and metal lines, will enable one to quantify the 
total Pop III He II emission and the Pop III spectral hardness (i.e. 
the ratio of He II to H ionizations) as a function of redshift. This 
will constrain the Pop III IMF and the evolution of the Pop III 
SFRD. Constraints on Pop III stars may also be improved with in¬ 


tensity mapping ofHg cooling lines which are produced during star 
formation Jciong et alj2013l) . Future work is required to determine 
how to best combine intensity mapping of different emission lines 
to put the tightest possible constraints on Pop III star formation. 

We estimated the Pop III He II 1640 A intensity mapping sig¬ 
nal from atomic cooling halos at z = 10 — 20 because atomic 
cooling halos are likely to be brighter and at lower redshift than 
minihalos. However, it may be possible to detect the signal from 
minihalos as well. Their signal at z ~ 20 could be nearly as high 
as the atomic cooling halo signal at z ~ 10 if the LW background 
were very low and/or ineffective in sterilizing minihalos (see Figure 
0}. We consider two examples of future intensity mapping observa¬ 
tions. In the first. Hell 1640 A is cross-correlated with CO(l-O) 
emission from galaxies and in the second He II is cross-correlated 
with 21 cm radiation from the IGM. Cross-correlation is necessary 
to remove contamination due to line emission from foreground and 
background galaxies. We find that in both examples the cross power 
spectrum can be measured with high signal-to-noise for Salpeter 
IMFs with mass ranges of (50 — 500)Mq and (1 — 500 )M 0 at 
z ~ 10. 

There are several caveats that could impact the effectiveness 
of this technique. First, there may be some degeneracy in the total 
He II signal and spectral hardness for Pop III metal-free stars with 
a particular IMF and stars which have a very low metallicity and a 
slightly different IMF (see fig. 5 in lSchaeren2003l) . There may also 
be degeneracy for different functional forms of the Pop III IMF. 
However, in any case, a very hard spectrum is likely to be strong 
evidence of a top-heavy Pop III IMF. Throughout we have assumed 
that the escape fraction of ionizing photons is much less than unity. 
We regard this as a reasonable assumption. If Hell-ionizing pho¬ 
tons were to leak out of galaxies, they would still be absorbed in the 
local IGM and produce only a small reduction in our signal (note 
that the Hell recombination time in the IGM at z > 10 is shorter 
than the Hubble time). Additionally, we stress that our estimates 
for the contribution from quasars are very uncertain. If the quasar 
signal is much higher than we estimate in section 6 , He II intensity 
mapping may be a better probe of high-redshift quasars than Pop 
III stars. This would still be very valuable; the inferred presence of 
a population of faint undetected miniquasars could be studied by 
correlations with the X-ray observations. Our estimate of the He II 
signal is also only intended to be a rough estimate. If the true signal 
is much lower than the values we compute, it may be necessary to 
design significantly more ambitious instruments than those we de¬ 
scribe for a high signal-to-noise detection. Despite all of the com¬ 
plications described above, because Pop III galaxies are so difficult 
to study individual l y, aside from observa tions of Supernovae (e.g. 
IWise & Abelll2005l ; fMesinger et alj|2006l) . He II intensity mapping 
may be the most direct method of constraining the Pop III IMF and 
the SFRD evolution at z > 10 in the foreseeable future. 
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